Establishing non-Abelian topological order in Gutzwiller projected Chern insulators 
via Entanglement Entropy and Modular 5-matrix 
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We use entanglement entropy signatures to establish non-Abelian topological order in a new class 
of ground states, the projected Chern-insulator wave functions. The simplest instance is obtained 
by Gutzwiller projecting a filled band with Chern number C = 2 which may also be viewed as 
the square of the band insulator Slater determinant. We demonstrate that this wave function is 
captured by the SU (2)2 Chern Simons theory coupled to fermions. In addition to the expected 
torus degeneracy and topological entanglement entropy, we also show that the modular iS-matrix, 
extracted from entanglement entropy calculations, provides direct access to the peculiar non-Abelian 
braiding statistics of Majorana fermions in this state. We also provide microscopic evidence for 
the field theoretic generalization, that the N th power of a Chern number C Slater determinant 
realizes the topological order of the SU(N)c Chern Simons theory coupled to fermions, by studying 
the SU(2)s and the SU(3)2 wave functions. An advantage of projected Chern insulator wave 
functions over lowest Landau level wave functions for the same phases is the relative ease with 
which physical properties, such as entanglement entropy, can be numerically calculated using Monte 
Carlo techniques. 

PACS numbers: 



It is known that quasiparticles may go beyond the 
conventional bosonic and fermionic statistics in two- 
dimensional many-body systems. Such systems represent 
a new class of phases not characterized by conventional 
order but by topological order—. One famous example is 
the Abelian quantum hall liquid for a fractionally filled 
Landau level in an interacting two-dimensional electron 
liquid 2 - and the corresponding Laughlin wave function^, 
where the quasiparticles exhibit fractional statistics. The 
low-energy effective theory is the U(l) Chern Simons 
(CS) theory and is relatively well understood 4 .. 

There has been much recent interest in an even more 
exotic generalization of statistics, non-Abelian statistics, 
displayed by phases with non-Abelian topological order. 
There is an increasing interest partially brought on by 
the recent preliminary evidence of Majorana fermions in 
superconductor-semiconductor junctions^ and their po- 
tentials for topological quantum memory and computa- 
tion realizations^—. Other examples include the v = 5/2 
fractional quantum hall effect^ and the Moore- Read Pfaf- 
fian stated. Although a general theory and classifica- 
tion of such states is still absent, there is a wide class 
effectively characterized by non-Abelian higher-level CS 
theories^ 1 -: the low-energy theory for filling v — k/N 
over k Landau levels is the SU(N)k CS theory, whose 
low energy effective theory is the SU (2) 2 CS theory. The 
elementary excitations are Ising anyons, whose braiding 
transforms the ground state, instead of just incur addi- 
tional phase factors as in the Abelian cases 1 ^. Its entan- 
glement spectrum has also been established as a powerful 
tool of identification 1 ^. For a universal topological quan- 
tum computation, the complexity of SU(2)s topological 
order with Fibonacci anyons (or above) is necessary 8 -, 
which may describe the fractional quantum hall plateau 
at v = 12/5. Examples of such phases are given by the 
Read-Rezayi states 1 ^ and generalizations to lattice spin 



liquid state s 15 ! 16 , however, their physical properties are 
relatively difficult to evaluate. 

For fractional quantum Hall liquids, extreme experi- 
mental conditions such as clean samples and large mag- 
netic fields are necessary to increase interaction and sup- 
press Landau level mixing. Yet it is later realized that 
the a Landau level structure is not necessary-^, as the 
quantum Hall effects is also present in Chern insulators 
- lattice band models without an external magnetic field 
but carry nontrivial reciprocal space Berry curvatures. It 
also offers a new arena for topological orders, and there 
have been various evidences for Abeliani 8 - - — and non- 
Abcliar>22r— fractional Chern insulators. 

Recall, the Laughlin wave functions^ ip ~ 
^[ (zi — Zj) m e~! 2fc ' / 4 , can be considered as the m th 
power of a integer Quantum Hall state of a particle 
with reduced charge, that fills the lowest Landau level. 
Previously, we confirmed that the lattice analog of 
this statement: the 771 th power of a Chern band wave 
function with unit Chern number %j) ~ ^c=l nas ^ ne 
topological order of a Laughlin state of order m 27 ' 28 . 
In this paper, we focus on the cases when the Chern 
number C > 1, which is unique to the lattice models and 
has no simple Landau level counterpart. Consistent with 
the field theory study in Ref. [29| and parton construction 
scheme proposed in Ref. H3, we suggest that the square 
(power N — 2) of C = 2 Chern band wave functions 
tp ~ V>c=2 are captured by the SU(2) 2 CS theory 
coupled to fermions and have the same quasiparticle 
braiding statistics as the Moore Read Pfaffian stated 

ip ~ pf (i-3i-)n0i- 2j)e- |zfc|2/4 . We verify such 
ground states are only three-fold degenerate. Especially, 
the wave-function diagnostic algorithm in Ref. l28| can be 
generalized to non-Abelian cases and particularly useful 
for many-body system where entanglement spectrum 
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is not available. With variational ansatz, physical 
measurable of these states are much simpler to calculate, 
thus we are able to extract the modular iS-matrix easily 
and determine the quantum dimension and quasiparticle 
statistics through topological entanglement entropy 
(TEE) 27 ' 28 ' 31 ' 32 , and prove the existence of non-Abelian 
quasiparticles. To our knowledge, this is the first direct 
numerical measurement of the modular 5-matrix and 
identification of a non-Abelian topological order wave 
function. We also extend our studies of ground-state 
degeneracy and entanglement to the generalized cases 
of C = 2,N = 3 (V> ~ ^0=2) and C* = 3,iV = 2 
(-0 ~ tpQ_ s ). These results imply the effective theory 
for the ./V th power of a C = k Chern insulator's band 
"0 ~ "0c=fe i s tne SU(N)k CS theory coupled to fermions, 
allowing non-Abelian statistics when N > 1 and k > 1. 
Besides, one may have access to the entire ground-state 
manifold with choices of boundary conditions of the 
parent Chern insulator. 

In the construction of a chiral spin liquid wave func- 
tion within the slave particle formalis m 27 ' 33 " —, the spin 
operator may be expressed as = /I/4., where ff and 
fl are annihilation operators for spin up and spin down 
fermions. If the / fermions fill up a band with Chern 
number C — 1 and is restricted to one fermion per site 
Hilbcrt space by Gutzwillcr projection, it will be a candi- 
date wave function for chiral spin liquid. Here we extend 
the situations to Chern number C > 1 to look for occur- 
rence of non-Abelian statistics. 

Chern Number C = 2 Model: To construct a two-band 
model with Chern number C = ±2, consider a tight- 
binding model on a two-dimensional square lattice with 
two orbitals on each lattice site labeled by / = 1, 2: 

H = E (-i) /+1 sW+E (^c] 2 c ll + h.c) 

<ij>,I <ij> 

+ & E (e 09 *4*i +*•<:■) (1) 

<<ik>> 

where Oij is the azimuthal angle for the vector connecting 
i and j. We rewrite the Hamiltonian in momentum space: 

H k = g(k) ■ 5 

gi = 2(cos k x — cos k y ) 

92 = 2<5 [cos(k x + k y ) — cos(fc x — k y )\ 

93 = 2(cos k x + cos k y ) (2) 

where Ui are Pauli matrices denoting the two orbitals. 
The model has a finite gap between the two bands with 
Chern number C = ±2 respectively. We have also veri- 
fied C = ±2 numerically by counting the number of chiral 
modes on the physical edges as well as within the entan- 
glement spectrum. Hereafter, we assume 6 = 1/V2 for a 
maximum gap to suppress the finite size effect. 

Torus Ground States of Gutzwiller Projected States 
and SU(2)2 CS Theory: Similar to the chiral spin liquid 
case2I, Eqn. [1] satisfies the Marshall sign rule, so 



the spin wave function may simply be express as the 
square of the fermionic wave function x 2 (z\ , z 2 • • • ) 
upto an unimportant phase, where x( z ij z 2 ■ ■ ■ ) is the 
ground-state wave function of lattice Chern insulator 
in Eqn. [T] and z = (r, I) contains both the position 
and orbital index. Its properties are the major focus 
of this paper. Note that x 2 ( z i> z i ' ' ' ) is 7r /2 rotational 
symmetric even though x(zi, z% ■ ■ • ) is not; although the 
7r/2 rotation symmetry is not essential to the topological 
properties, it will be convenient for determining them. 

Physical quantities related with the wave func- 
tions may be calculated with variational Monte Carlo 
method 3 ^. We study the wave functions on a torus with 
periodic or antiperiodic boundary conditions along the x 
and y directions in Eqn. [TJ Note that the squared wave 
function \ 2 { z ii z 2 • • • ) is insensitive to such boundary 
conditions, hence all four such wave functions are equally 
good candidate ground states, denoted as I'&^t&j, > where 
&x.y = 0, it. To understand their linear dependence and 
the ground-state manifold, we calculate the overlaps be- 
tween them with variational Monte Carlo method on a 
12 x 12 system: 

(OOItttt) = a 

(OttIttO) = a 

(0tt|00) = (tt0|00) = /3 

(OttItttt) = (ttOItttt) = P' (3) 

Numerically we find to very high accuracy that a = 
a' = -0.086 and (} = (}' = 0.457. We may construct a 
"generalized" projection operator: 




|tt0) \ t /1 a P p \ ( (vr0| \ 

[Ott) a 1 p p \ (Ovrj 

\tttt) P p 1 a \ (tttt\ 

|00) / V p p a 1 / V (00| / 



(4) 

which projects to the ground state manifold. Due to the 
non-orthogonality between the basis states, the eigenval- 
ues of P contains one 0, so the corresponding eigenstate 
is projected out: 

P[\tttt > -| Ott > — |tt0 > +|00 >] = (5) 

where we have used 2/3 = 1 + a (true to high numerical 
accuracy). Thus the ground states are only three- fold 
degenerate on a torus: 

\F x = l,F y = l> ~ (|00 > +|0tt > +|tt0 > +\tttt >) 
\F X = l,F y = -l> ~ (|00 > -| Ott > +|tt0 > -\tttt >) 
\F X =-l,F y = l> ~ (|00 > +|0tt > -|tt0 > -\ttit >) 

(6) 

upto phase and normalization, where we have used the 
basis of flux threading operators F x and F y to label these 



3 



states. The ground-state degeneracy and linear depen- 
dence is consistent with the SU(2) 2 CS theory, as shown 
in Rcf. 41. In addition to these simple choices of bound- 
ary conditions, there are other possible constructions for 
ground states, where the wave function is a product of 
two fermionic wave functions with different boundary 
conditions. These are shown to be linearly dependent 
on the wave functions above^i. 

Topological Entanglement Entropy: To obtain further 
information on the wave functions, we use entanglement 
as a tool and extract their TEE on a 6 x 6 system with a 
method by Kitaev and Preskill^I, which effectively can- 
cels out the boundary and corner contributions and ex- 
poses the topological term when the typical length scale 
is much longer than the correlation length. Although the 
smallest length scale is only 2 lattice spacings for the sys- 
tem size we study, it is still longer than the correlation 
length £ ~ 0.5 lattice spacing. In addition, the corre- 
sponding wave-function overlaps suggest that the residue 
of Eqn. [5] is just ~ 1.3%, thus the three-fold degeneracy 
in Eqn. [5] is still a good approximation. These facts sug- 
gest that the system size is large enough to suppress finite 
size effect and maintain topological properties. Through- 
out we focus on the Renyi entropy S2 for algorithmic 
simplification s 27 ' 28 ' 37 . 

We find that the TEE of \&x$y > for a topologically 
trivial disc-shape cut is 7 = 0.85 ± 0.08, in fairly good 
consistency with the theoretical value Xsu<2) 2 = 1°§2 ~ 
0.693 for the SU{2) 2 CS theory. Note that for an Abelian 
topological order with D 2 = 3 degenerate ground states 
on a torus, the expected TEE 7 = logD ~ 0.549 is fur- 
ther away from the calculated value and unlikely to de- 
scribe the wave functions. It implies non-Abelian statis- 
tics. 

Modular S-Matrix from Entanglement Entropy: Next 
we extract the quasiparticlc braiding information using 
entanglement. Following the algorithm in Ref. [28|, we (i) 
calculate the TEE 7' for partitioning the torus into two 
cylinders along the y direction for linear combinations 
of wave functions, see Fig. [TJ inset, (ii) search for the 
states with minimum entanglement entropy (maximum 
TEE 7') and identify them as the quasiparticle Wilson 
loop states, and (iii) study the 7r/2 rotation transforma- 
tion on them. The resulting TEE 27 - 7' for |$i >= 
cos0|O7r > +sin(9|7rO >, |$ 2 >= sin6»|00 > -cos6>|7rO >, 
and: 



|$ 3 > = (sin(9 + 0.7915cos(9)|00 > 

- (sin 6> + 0.4697 cos 6) |tt0 > 

- 1.2623 cos 9 1 Ott > 



(7) 



for selected values of 9 is shown in Fig. [TJ Fig. [3^ and 
Fig. [2Jd, respectively. 

Then we have obtained approximately the three or- 
thogonal quasiparticlc Wilson loop states corresponding 
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FIG. 1: (Inset) Kitaev Preskill scheme for extracting TEE 
by partition the torus into subsystems A, B and C and con- 
structing A = -Sabc + Sab + Sbc + Sac -Sa — Sb — Sc = 
—2Sa + "2Sab — Sabct—- Note here regions C and AB en- 
circle the torus leading to TEE of 7', which is ground state 
dependent, unlike the TEE for region A (or B) which is al- 
ways 7 28 ' 37 . The resulting A = 27 — 7' is plotted for the 
linear combinations of wave functions |$i >. The solid curve 
is expected value for the SU (2)2 CS theory. 



to the minima of 27 — 7—: 

|Hi > = -|tt0 > -|00 > 
IS2 > = — |tt0 > +|00 > 

|H 3 > = 0.7915|00 > -0.4697|tt0 > - 1.2623 |0?r >(8) 

For a 7r/2 rotation symmetric system, the modular S- 
matrix's action in the |Sj- > basis is a 7r/2 rotation. Note 
the transformations |-7r0 ><-> |07r > and |00 |00 > 
under 7r/2 rotation, we may transform the basis to Eqn. 
[S]with the help of Eqn. [3]to properly treat orthogonality: 



0.627 0.610 0.484 
0.610 0.000 -0.792 
0.484 -0.792 0.372 



(9) 



The rough estimation of the modular 5-matrix in Eqn. 
[5]already suggests important qualitative results for quasi- 
particle braiding. Especially, the phase for the Jj 'th entry 
of the (S-matrix is the phase factor during the process 
of creating the i'th and j'th quasiparticlc, braiding the 
quasiparticles before annihilating them with each anti- 
quasiparticles and returning to the ground states. While 
the quasiparticles corresponding to the first and third 
rows and columns obey Abelian bosonic or fermionic 
statistics for a trivial phase upon braiding, the zero value 
entry for the second quasiparticle is a signature to its non- 
Abelian self statistics. Indeed, for the Majorana fcrmion 
in the SU(2)2 CS theory, one fcrmion composes a pair 
of Majorana fermions c = 71 + 272, however, when one 
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FIG. 2: TEE A = 27 - 7' for linear combinations of wave 
functions: (a) |$2 >, and (b) $3 >. The solid curves are 
expected values for the 5*17(2)2 CS theory. Note, the presence 
of only two minima with A ~ indicates that two of the three 
quasiparticles are Abelian, while the third is non-Abelian. 



of the Majorana fcrmions braids with another Majorana 
fermion and picks an additional n phase, it changes the 
fermion annihilation operator to creation operator and 
vice versa, and fails to return to the ground state after- 
wards, thus the corresponding zero entry in the modular 
iS-matrix. 

As a comparison, the ideal iS-matrix for the SU(2)2 
CS theory is: 

/ 0.5 0.707 0.5 \ 
S = 0.707 -0.707 (10) 
\ 0.5 -0.707 0.5 / 

roughly consistent with Eqn. |H1 Note that although the 
|Sj > states in Eqn. [5] are estimates, the qualitative fea- 
tures from the corresponding 5-matrix are more robust. 

For a more detailed comparison, we note the SU(2)2 
CS theory predicts the following connection between the 
ground state wave functions of Eqn|5] and the anyons of 



this phased: 

\F X = l,Fy = l> = (\ly>+\i> y >)/V2 
\F X = \,Fy = -l> = {\ly> -|V>„ >)/V2 

\F X =-l,Fy = l> = \Uy > (11) 

where \l y >, \ip y > and \<r y > are the Wilson loop states 
along the y direction for the identity, fermionic and non- 
Abelian Ising anyonic quasiparticles, respectively. Then 
for an arbitrary ground state: 

27-V = -bgfc>2/dA (12) 

where dj is the individual quantum dimension and pj 
is the statistical weight in the wave function for the j 
quasiparticles^. Combining Eqn. [3J Eqn. [5] and Eqn. 
1111 with relative phase between wave functions therein, 
one may expand the l^z^y > states in the \l y >, \ip y >, 
\<j y > basis. The corresponding values of 27 — 7' (6) given 
by Eqn. [T^] and d\ = = l,d a = y2 are shown in 
Fig. [1] and Figj2] as the solid curves, and fit well with 
the numerical results. These individual quantum dimen- 
sions dj are also implicit in the values of the minima as 
27 — 7^ = log(<i|), which directly follows from Eqn. [T2] 
Especially, d a = a/2 suggests the o quasiparticle obeys 
non-Abelian statistics. 

Other non-Abelian States (i) SU(3)2'- Such construc- 
tion of wave functions may be generalized to more com- 
plicated non-Abelian cases. As another example, for the 
wave functions x 3 (zi > z 2 • • • ) , we can construct nine can- 
didate ground-state wave functions | c E> a ; ( l > ; / > with simple 
boundary conditions $ x .y = 0,±27r/3. Repeating the 
procedure for Eqn [3] and Eqn. [4j we obtain its "gener- 
alized" projection operator P' on a 12 x 12 torus, see 
Ref. 41. Indeed, P' only has six non-zero eigenvalues, 
therefore the ground states are six-fold degenerate. The 
six corresponding eigenstates' 7r/2 rotation eigenvalues 
are ±1, ±1 and ±i. These results arc consistent with 
the SU(3) 2 CS theory. In addition, the TEE for a disc- 
shaped cut on a 6 x 6 torus is 7 ~ 1.264 ± 0.073, con- 
sistent with the theoretical value of D = \Js(5 + \fb)/2 
and 7s(7(3) 2 = l°g D ~ 1.19. While we consider the above 
evidences sufficient, we leave further verifications as TEE 
ground-state dependence and linear dependence of more 
complex boundary constructions to future works. 

(ii) SU(2)z in close connection to the Read-Rezayi 
state: In the same spirit, to construct wave functions 
for the Fibonacci anyons characterized by the SU(2)s 
CS theory, we may construct the following tight-binding 
model on a triangular lattice with two orbitals on each 
site labeled by I = 1,2 in analogy to Eqn. [TJ 

H = ]T {-l) I+1 c) lCu+ Y. (e i3( S 2Cil + />.c.) 

<ij>,I <ij> 

+ S (e i3eik cl 2 c a + h.c}j (13) 

<<ifc>> 
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where 9ij is the azimuthal angle for the vector connecting 
i and j. The model is a two-band Chern insulator with 
Chcrn number C = ±3 for 5^0. Similar construction 
may have potentials for the construction of even higher 
Chcrn number bands, and a systematic scheme to pro- 
duce arbitrary Chern number bands has been studied in 
Ref.lH We choose nine different boundary conditions in 
Eqn. [T3]for the construction of X 2 ( z ii z i ' ' ' ) on a 12 x 12 
torus. Our results show that only four of the nine eigen- 
values of the corresponding "generalized" projection op- 
erator have significantly finite values, so the ground state 
is only four-fold degeneracy, consistent with the SU (2)s 
CS theory. 

These examples indicate that the iV th power of a Chern 
insulator with Chern number C = k is characterized by 
the SU(N)k CS theory (coupled to fcrmions), thus pro- 



viding a novel route to construct wave functions for non- 
Abelian topological order, and the power of entanglement 
in the diagnosis and study of these wave functions. 

In conclusion we have introduced a class of lattice wave 
functions for non-Abelian topological phases that (i) are 
readily generalized to capture SU(N)k topological order 
and (ii) easily generate all degenerate ground states and 
(iii) whose properties can be compute with using monte- 
carlo techniques. Further the presence of such natural 
lattice wave functions holds promise that such states may 
be realized in the context of fractional Chcrn insulators. 
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Strictly speaking this wave function is not the same topo- 
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by the pure SU (2)2 CS theory; although it shares the 
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braiding encoded in the modular S-matrix with the Pfaf- 
fian state, it differs at the level of the quasiparticle spin 
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See Ref. [3j| for detailed discussion. We do not focus on 
such distinctions because these states behave identical for 
the properties we study. 
See supplementary material online. 
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Appendix A: Parton boundary conditions for Chern 
Simons Theories 



The SU(N) k CS theory is : 
£su(N) k = -^£p.v\Tr (^a^dvax + -a^a^a^j , a M <G SU(N) 

In this appendix we first study the boundary condi- 
tions and linear dependence for various wave functions 
through parton construction for SV "(2)2 CS theory Af- 
ter choosing gauge ao = and a = a ■ Diag(l, —1), we 
have: 

C-su{2) 2 = -^ a ^d v a\ 

denoting low energy degrees of freedom u = a 1 L 1 /27r and 
v = a2l/2/27r, both defined modulo 1, the effective action 
becomes: 

C e ff = 8ttvu 
therefore the canonical commutation relation: 



which have the same boundary condition for both parton 
wave functions x, and: 



\u,v 



8" 



so the conjugate momentum of u is 8irv, with (u,v) ~ 
(—it, —v) due to global SU(2) gauge transformations. 
Then if the wave function in the v coordinate is: 

4>(v) = c„<5(4t; — n) 

which is essentially the Wilson loop states in the y di- 
rection. Similarly, the wave function in the u coordinate 
is: 



4>(u) = y^c„e' 



As a result, in the spirit of parton construction, we 
may write down the following wave functions: 



\h > = 



and: 



0(0,0) + 0(1/2,0) + 720(1/4,0)] /2 
0(0, 1/2) + 0(1/2, 1/2) - \/20(l/4, 1/2)] /2 
[0(0, 1/4) -0(1/2, 1/4)} /V2 



(|1„>-|^>)/V2 = [0(1/4,0) -0(1/4, 1/2)] fy/2 
{\l v >+\^ y >)/V2 = [0(0,0) + 0(1/2,0) 

+ 0(0, 1/2) + 0(l/2,l/2)]/2 

here eight wave functions with different parton boundary 
conditions are involved: 



0(0,0) 
0(0,1/2) 
0(1/2,0) 

0(0,0) 



2 

Xo,o 

Xtt.O 



0(0,1/4) 
(1/2,1/4) 

0(1/4,0) 
(1/4,1/2) 



X0,7r/2X0,-7r/2 
Xtt.it /2Xtt,-tt /2 
Xtt/2,0X-tt/2,0 
Xtt /2,ttX— tt/2,tt 



where the boundary conditions for two partons are dif- 
ferent. 

However, these wave functions are not completely lin- 
ear independent. Comparing the rotation eigenstates and 
the Wilson loop states, we obtain the following linear de- 
pendence: 

[0(0, 1/4) - 0(1/4, 0) + 0(1/4, 1/2) - 0(1/2, 1/4)] /2 
= [0(0, 1/2)- 0(1/2, 0)]/V5 

[0(0, 1/4) + 0(1/4, 0) + 0(1/4, 1/2) + 0(1/2, 1/4)] /2 
= [0(0, 0) + 0(1/2, 0) + 0(1/2, 0) + 0(1/2, 1/2)] /2 
[0(0, 1/4) + 0(1/4, 0) - 0(1/4, 1/2) - 0(1/2, 1/4)] /2 
= [0(0,0) -0(1/2, 1/2)] /V2 

0(0, 1/4) - 0(1/4, 0) - 0(1/4, 1/2) + 0(1/2, 1/4) = 



So the latter four wave functions may be fully established 
from the former four. In addition: 

0(0, 0) - 0(0, 1/2) - 0(1/2, 0) + 0(1/2, 1/2) = 

which we have numerically verified as Eqn. [5] in the main 
text of the paper. To verify the other four linear de- 
pendence relations, we obtain numerically the following 
overlaps on a 12 x 12 torus between the wave functions 
with relevant boundary conditions in Eqn. [TJ 

(00|0cr) = (00|ct0) ~ (tttt\tt(t) = (ttttIcttt) 
~ (ttOIttct) ~ (07r|0cr) ~ 0.8403 

(00|7rcr) = (OOIcttt) ~ (ttttIOct) = (ittt\<t0) 
~ (OttIttct) ~ (7r0|0cr) ~ 0.0726 

where | Oct) , |7rcr), |cr0) and \air) correspond to 0(0,1/4), 
0(1/2,1/4), 0(0,1/2) and 0(1/4,1/2), respectively. 
Combining these results with Eqn. [31 one can easily re- 
produce the rest four relations. 

Overall, there are five constraints so leave only three 
linearly independent wave functions, consistent with the 
three-fold ground-state degeneracy in a SU (2)2 CS the- 
ory. 

The above discussion may be generalized to parton 
constructions of other SU(N)k CS theories. It indicates 
certain relations of wave functions constructed from vari- 
ous boundary conditions, and implies a method to verify 
or identify the underlying topological theory of a set of 
parton wave functions just by calculating their overlaps. 
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As a second example, we calculate the wave-function 
overlap for x 3 ( z i: z 2 • • • ) of Eqn. U on a 12 x 12 sys- 
tem. In the basis of (|00) , |o4f), \0=^), |^0>, \^^-), 

llf^r)> l^ )' I^Tlf)> I ^^r))- the correspond- 
ing "generalized" projection operator is: 



/I 

a 



P' 



a 
a 
1 



b 

ac* 
ac- 



he be* 

ac be* 

be ac* 
be* be 
b ac* be be* ac 

\ b be* ac be b 



a 

be* 
be 

1 ac ac* 
ac* 1 ac 
ac ac* 1 
a be* be 
b 

_ * 
ac 



b a b b \ 

be* be ac be 

be be* be* ac* 

a be be* 

be ac* b 

be* b ac 

1 ac* ac 

ac 1 ac* 

ac* ac 1 ) 



where a ~ 0.4213, b ~ 0.0768 ~ 0.5 - a, c ~ e l27r / 3 from 
numerical calculations. P' only has six non-zero eigenval- 
ues, therefore the ground states are six-fold degenerate. 
The corresponding cigenstates' 7r/2 rotation eigenvalues 
are ±1, ±1 and ±i, consistent with the SU (3)2 CS theory. 
Especially, by comparing the rotation eigenvalues of the 
modular 5-matrix with those of P' eigenstates, one may 
obtain useful clues for the relations of the quasiparticlc 
Wilson loop states with l^x^y >■ 



